Evaluating differences in milk production, reproductive performance, and survival associated with vaginal discharge characteristics and fever in postpartum dairy cows.

The objective was to assess differences in productive and reproductive performance, and survival associated with vaginal discharge characteristics and fever in postpartum dairy cows located in Western and Southern states of the U.S.A. This retrospective cohort study included data from 3 experiments conducted in 9 dairies. Vaginal discharge was evaluated twice within 12 DIM and scored on a 5-point scale. The highest score observed for each cow was used for group assignment (VD group) as follows: VD 1 and 2 (VD 1/2; n = 1,174) = clear mucus/lochia with or without flecks of pus; VD 3 (n = 1,802) = mucopurulent with < 50% pus; VD 4 (n = 1,643) = mucopurulent with ≥50% of pus or non-fetid reddish/brownish mucous, n = 1,643; VD 5 = fetid, watery, and reddish/brownish, n = 1,800. All VD 5 cows received treatment according to each herd's protocol. Rectal temperature was assessed in a subset of VD 5 cows, and subsequently divided into Fever (rectal temperature ≥39.5°C; n = 334) and NoFever (n = 558) groups. A smaller proportion of cows with VD 5 (67.6%) resumed ovarian cyclicity compared with VD 1/2 (76.2%) and VD 4 (72.9%) cows; however, a similar proportion of VD5 and VD 3 (72.6%) cows resumed ovarian cyclicity. A smaller proportion of VD 5 (85.8%) cows received at least one artificial insemination (AI) compared with VD 1/2 (91.5%), VD 3 (91.0%), or VD 4 (91.6%) cows. Although we did not detect differences in pregnancy at first AI according to VD, fewer cows with VD 5 (64.4%) were pregnant at 300 DIM than cows with VD 1/2 (76.5%), VD 3 (76.2%), or VD 4 (74.7%). Hazard of pregnancy by 300 DIM was smaller for VD 5 compared with VD 1/2, VD 3, or VD 4 cows. A greater proportion of VD 5 cows were removed from the herd within 300 DIM compared with other VD groups. There was 760 kg lesser milk production within 300 DIM for VD 5 compared with VD 2, VD 3, and VD 4, whereas VD 2, VD 3, and VD 4 had similar milk production. We did not detect an association between fever at diagnosis of VD 5 and reproductive performance or milk production. A greater proportion of VD 5 cows without fever were removed from the herd by 300 DIM compared with VD 5 cows with fever. Differences in productive and reproductive performance, and removal of the herd were restricted to fetid, watery, and reddish/brownish vaginal discharge, which was independent of fever.

Autonomic cardiac function in children and adolescents with long COVID: a case-controlled study.

Although the mechanisms underlying the pathophysiology of long COVID condition are still debated, there is growing evidence that autonomic dysfunction may play a role in the long-term complications or persisting symptoms observed in a significant proportion of patients after SARS-CoV-2 infection. However, studies focused on autonomic dysfunction have primarily been conducted in adults, while autonomic function has not yet been investigated in pediatric subjects. In this study, for the first time, we assessed whether pediatric patients with long COVID present abnormalities in autonomic cardiac function. Fifty-six long COVID pediatric patients (mean age 10.3 ± 3.8 y) and 27 age-, sex-, and body surface area-matched healthy controls (mean age 10.4 ± 4.5y) underwent a standard 12-lead electrocardiography (ECG) and 24-h ECG Holter monitoring. Autonomic cardiac function was assessed by time-domain and frequency-domain heart rate variability parameters. A comprehensive echocardiographic study was also obtained by two-dimensional echocardiography and tissue Doppler imaging. Data analysis showed that pediatric patients with long COVID had significant changes in HRV variables compared to healthy controls: significantly lower r-MSSD (root mean square of successive RR interval differences, 47.4 ± 16.9 versus 60.4 ± 29.1, p = 0.02), significant higher values VLF (very low frequency, 2077.8 ± 1023.3 versus 494.3 ± 1015.5 ms, p = 0.000), LF (low frequency, 1340.3 ± 635.6 versus 354.6 ± 816.8 ms, p = 0.000), and HF (high frequency, 895.7 ± 575.8 versus 278.9 ± 616.7 ms, p = 0.000). No significant differences were observed between the two groups both in systolic and diastolic parameters by echocardiography.  Conclusion: These findings suggest that pediatric patients with long COVID have an imbalance of cardiac autonomic function toward a relative predominance of parasympathetic tone, as already reported in adult patients with long COVID. Further studies are needed to clarify the clinical significance of this autonomic dysfunction and demonstrate its role as a pathophysiological mechanism of long COVID, paving the way for effective therapeutic and preventive strategies. What is Known: • Long Covid in children has been described globally, but studies have mostly focused on collecting the temporal evolution of persisting symptoms. What is New: • Cardiac autonomic imbalance toward a relative predominance of parasympathetic tone is a mechanism underlying Long Covid in children, as also described in adults.

A systematic review of influenza virus in water environments across human, poultry, and wild bird habitats.

Influenza, a highly contagious acute respiratory disease, remains a major global health concern. This study aimed to comprehensively assess the prevalence of influenza virus in different aquatic environments. Using 43 articles from four databases, we thoroughly examined water matrices from wastewater treatment plants (WTPs) and other human environments, as well as poultry habitats and areas frequented by migratory wild birds. In WTP influents (10 studies), positivity rates for influenza A ranged from 0.0 % to 97.6 %. For influenza B (8 studies), most studies reported no positivity, except for three studies reporting detection in 0.8 %, 5.6 %, and 46.9 % of samples. Within poultry habitats (13 studies), the prevalence of influenza A ranged from 4.3 % to 76.4 %, while in environments frequented by migratory wild birds (11 studies), it ranged from 0.4 % to 69.8 %. Geographically, the studies were distributed as follows: 39.5 % from the Americas, 18.6 % from Europe, 2.3 % from South-East Asia and 39.5 % from the Western Pacific. Several influenza A subtypes were found in water matrices, including avian influenza (H3N6, H3N8, H4N1, H4N2, H4N6, H4N8, H5N1, H5N8, H6N2, H6N6, H7N9, H0N8, and H11N9) and seasonal human influenza (H1N1 and H3N2). The existing literature indicates a crucial requirement for more extensive future research on this topic. Specifically, it emphasizes the need for method harmonization and delves into areas deserving of in-depth research, such as water matrices pertaining to pig farming and prevalence studies in low-income countries.

Genome-enabled prediction of indicator traits of resistance to gastrointestinal nematodes in sheep using parametric models and artificial neural networks.

This study aimed to assess the predictive ability of parametric models and artificial neural network method for genomic prediction of the following indicator traits of resistance to gastrointestinal nematodes in Santa Inês sheep: packed cell volume (PCV), fecal egg count (FEC), and Famacha© method (FAM). After quality control, the number of genotyped animals was 551 (PCV), 548 (FEC), and 565 (FAM), and 41,676 SNP. The average prediction accuracy (ACC) calculated by Pearson correlation between observed and predicted values and mean squared errors (MSE) were obtained using genomic best unbiased linear predictor (GBLUP), BayesA, BayesB, Bayesian least absolute shrinkage and selection operator (BLASSO), and Bayesian regularized artificial neural network (three and four hidden neurons, BRANN_3 and BRANN_4, respectively) in a 5-fold cross-validation technique. The average ACC varied from moderate to high according to the trait and models, ranging between 0.418 and 0.546 (PCV), between 0.646 and 0.793 (FEC), and between 0.414 and 0.519 (FAM). Parametric models presented nearly the same ACC and MSE for the studied traits and provided better accuracies than BRANN. The GBLUP, BayesA, BayesB and BLASSO models provided better accuracies than the BRANN_3 method, increasing by around 23% for PCV, and 18.5% for FEC. In conclusion, parametric models are suitable for genome-enabled prediction of indicator traits of resistance to gastrointestinal nematodes in sheep. Due to the small differences in accuracy found between them, the use of the GBLUP model is recommended due to its lower computational costs.

Detection and full genomic sequencing of rare hepatitis E virus genotype 4d in Italian wastewater, undetected by clinical surveillance.

Hepatitis E is a liver disease caused by the hepatitis E virus (HEV), primarily transmitted through contaminated water or food. There are four different HEV genotypes in humans, with genotypes 1 and 2 being the most widespread. Genotypes 3 and 4 are found in animals and can also infect humans. Genotype 4 is prevalent in Asia, mainly in China. In Italy, only one outbreak of HEV-4 has been documented, which occurred in 2011, involving five patients. In 2013, HEV G4 was also detected in a pig farm. Since then, no further evidence of HEV genotype 4 has been found in the country. This study describes the first detection of HEV genotype 4, subtype d, in wastewater in central Italy, despite a lack of any clinical case reported in the area. By using a multiplex PCR protocol and two sequencing strategies, Illumina and ONT, the virus's complete genome was sequenced and characterized as subtype 4d. These findings shed light on the potential of environmental surveillance for infectious agents to improve our understanding of epidemiology and support public health efforts.

Single-cell DNA methylome and 3D multi-omic atlas of the adult mouse brain.

Cytosine DNA methylation is essential in brain development and is implicated in various neurological disorders. Understanding DNA methylation diversity across the entire brain in a spatial context is fundamental for a complete molecular atlas of brain cell types and their gene regulatory landscapes. Here we used single-nucleus methylome sequencing (snmC-seq3) and multi-omic sequencing (snm3C-seq)1 technologies to generate 301,626 methylomes and 176,003 chromatin conformation-methylome joint profiles from 117 dissected regions throughout the adult mouse brain. Using iterative clustering and integrating with companion whole-brain transcriptome and chromatin accessibility datasets, we constructed a methylation-based cell taxonomy with 4,673 cell groups and 274 cross-modality-annotated subclasses. We identified 2.6 million differentially methylated regions across the genome that represent potential gene regulation elements. Notably, we observed spatial cytosine methylation patterns on both genes and regulatory elements in cell types within and across brain regions. Brain-wide spatial transcriptomics data validated the association of spatial epigenetic diversity with transcription and improved the anatomical mapping of our epigenetic datasets. Furthermore, chromatin conformation diversities occurred in important neuronal genes and were highly associated with DNA methylation and transcription changes. Brain-wide cell-type comparisons enabled the construction of regulatory networks that incorporate transcription factors, regulatory elements and their potential downstream gene targets. Finally, intragenic DNA methylation and chromatin conformation patterns predicted alternative gene isoform expression observed in a whole-brain SMART-seq2 dataset. Our study establishes a brain-wide, single-cell DNA methylome and 3D multi-omic atlas and provides a valuable resource for comprehending the cellular-spatial and regulatory genome diversity of the mouse brain.

Brain-wide correspondence of neuronal epigenomics and distant projections.

Single-cell analyses parse the brain's billions of neurons into thousands of 'cell-type' clusters residing in different brain structures1. Many cell types mediate their functions through targeted long-distance projections allowing interactions between specific cell types. Here we used epi-retro-seq2 to link single-cell epigenomes and cell types to long-distance projections for 33,034 neurons dissected from 32 different regions projecting to 24 different targets (225 source-to-target combinations) across the whole mouse brain. We highlight uses of these data for interrogating principles relating projection types to transcriptomics and epigenomics, and for addressing hypotheses about cell types and connections related to genetics. We provide an overall synthesis with 926 statistical comparisons of discriminability of neurons projecting to each target for every source. We integrate this dataset into the larger BRAIN Initiative Cell Census Network atlas, composed of millions of neurons, to link projection cell types to consensus clusters. Integration with spatial transcriptomics further assigns projection-enriched clusters to smaller source regions than the original dissections. We exemplify this by presenting in-depth analyses of projection neurons from the hypothalamus, thalamus, hindbrain, amygdala and midbrain to provide insights into properties of those cell types, including differentially expressed genes, their associated cis-regulatory elements and transcription-factor-binding motifs, and neurotransmitter use.

Conserved and divergent gene regulatory programs of the mammalian neocortex.

Divergence of cis-regulatory elements drives species-specific traits1, but how this manifests in the evolution of the neocortex at the molecular and cellular level remains unclear. Here we investigated the gene regulatory programs in the primary motor cortex of human, macaque, marmoset and mouse using single-cell multiomics assays, generating gene expression, chromatin accessibility, DNA methylome and chromosomal conformation profiles from a total of over 200,000 cells. From these data, we show evidence that divergence of transcription factor expression corresponds to species-specific epigenome landscapes. We find that conserved and divergent gene regulatory features are reflected in the evolution of the three-dimensional genome. Transposable elements contribute to nearly 80% of the human-specific candidate cis-regulatory elements in cortical cells. Through machine learning, we develop sequence-based predictors of candidate cis-regulatory elements in different species and demonstrate that the genomic regulatory syntax is highly preserved from rodents to primates. Finally, we show that epigenetic conservation combined with sequence similarity helps to uncover functional cis-regulatory elements and enhances our ability to interpret genetic variants contributing to neurological disease and traits.

Hotspots of Somatic Genetic Variation in Pituitary Neuroendocrine Tumors.

The most common genetic drivers of pituitary neuroendocrine tumors (PitNETs) lie within mutational hotspots, which are genomic regions where variants tend to cluster. Some of these hotspot defects are unique to PitNETs, while others are associated with additional neoplasms. Hotspot variants in GNAS and USP8 are the most common genetic causes of acromegaly and Cushing's disease, respectively. Although it has been proposed that these genetic defects could define specific clinical phenotypes, results are highly variable among studies. In contrast, DICER1 hotspot variants are associated with a familial syndrome of cancer predisposition, and only exceptionally occur as somatic changes. A small number of non-USP8-driven corticotropinomas are due to somatic hotspot variants in USP48 or BRAF; the latter is a well-known mutational hotspot in cancer. Finally, somatic variants affecting a hotspot in SF3B1 have been associated with multiple cancers and, more recently, with prolactinomas. Since the associations of BRAF, USP48, and SF3B1 hotspot variants with PitNETs are very recent, their effects on clinical phenotypes are still unknown. Further research is required to fully define the role of these genetic defects as disease biomarkers and therapeutic targets.

Clinical Spectrum of USP8 Pathogenic Variants in Cushing's Disease.

Cushing's disease (CD) is a life-threatening condition with a challenging diagnostic process and scarce treatment options. CD is caused by usually benign adrenocorticotrophic hormone (ACTH)-secreting pituitary neuroendocrine tumors (PitNETs), known as corticotropinomas. These tumors are predominantly of sporadic origin, and usually derive from the monoclonal expansion of a mutated cell. Somatic activating variants located within a hotspot of the USP8 gene are present in 11-62% of corticotropinomas, making USP8 the most frequent genetic driver of corticotroph neoplasia. In contrast, other somatic defects such as those affecting the glucocorticoid receptor gene (NR3C1), the BRAF oncogene, the deubiquitinase-encoding gene USP48, and TP53 are infrequent. Moreover, patients with familial tumor syndromes, such as multiple endocrine neoplasia, familial isolated pituitary adenoma, and DICER1 rarely develop corticotropinomas. One of the main molecular alterations in USP8-driven tumors is an overactivation of the epidermal growth factor receptor (EGFR) signaling pathway, which induces ACTH production. Hotspot USP8 variants lead to persistent EGFR overexpression, thereby perpetuating the hyper-synthesis of ACTH. More importantly, they condition a characteristic transcriptomic signature that might be useful for the clinical prognosis of patients with CD. Nevertheless, the clinical phenotype associated with USP8 variants is less well defined. Hereby we discuss the current knowledge on the molecular pathogenesis and clinical picture associated with USP8 hotspot variants. We focus on the potential significance of the USP8 mutational status for the design of tailored clinical strategies in CD.

Differential effects of follicle-stimulating hormone glycoforms on the transcriptome profile of cultured rat granulosa cells as disclosed by RNA-seq.

It has been documented that variations in glycosylation on glycoprotein hormones, confer distinctly different biological features to the corresponding glycoforms when multiple in vitro biochemical readings are analyzed. We here applied next generation RNA sequencing to explore changes in the transcriptome of rat granulosa cells exposed for 0, 6, and 12 h to 100 ng/ml of four highly purified follicle-stimulating hormone (FSH) glycoforms, each exhibiting different glycosylation patterns: human pituitary FSH18/21 and equine FSH (eqFSH) (hypo-glycosylated), and human FSH24 and chinese-hamster ovary cell-derived human recombinant FSH (recFSH) (fully-glycosylated). Total RNA from triplicate incubations was prepared from FSH glycoform-exposed cultured granulosa cells obtained from DES-pretreated immature female rats, and RNA libraries were sequenced in a HighSeq 2500 sequencer (2 × 125 bp paired-end format, 10-15 × 106 reads/sample). The computational workflow focused on investigating differences among the four FSH glycoforms at three levels: gene expression, enriched biological processes, and perturbed pathways. Among the top 200 differentially expressed genes, only 4 (0.6%) were shared by all 4 glycoforms at 6 h, whereas 118 genes (40%) were shared at 12 h. Follicle-stimulating hormone glycocoforms stimulated different patterns of exclusive and associated up regulated biological processes in a glycoform and time-dependent fashion with more shared biological processes after 12 h of exposure and fewer treatment-specific ones, except for recFSH, which exhibited stronger responses with more specifically associated processes at this time. Similar results were found for down-regulated processes, with a greater number of processes at 6 h or 12 h, depending on the particular glycoform. In general, there were fewer downregulated than upregulated processes at both 6 h and 12 h, with FSH18/21 exhibiting the largest number of down-regulated associated processes at 6 h while eqFSH exhibited the greatest number at 12 h. Signaling cascades, largely linked to cAMP-PKA, MAPK, and PI3/AKT pathways were detected as differentially activated by the glycoforms, with each glycoform exhibiting its own molecular signature. These data extend previous observations demonstrating glycosylation-dependent differential regulation of gene expression and intracellular signaling pathways triggered by FSH in granulosa cells. The results also suggest the importance of individual FSH glycoform glycosylation for the conformation of the ligand-receptor complex and induced signalling pathways.

Single-cell DNA methylation and 3D genome architecture in the human brain.

Delineating the gene-regulatory programs underlying complex cell types is fundamental for understanding brain function in health and disease. Here, we comprehensively examined human brain cell epigenomes by probing DNA methylation and chromatin conformation at single-cell resolution in 517 thousand cells (399 thousand neurons and 118 thousand non-neurons) from 46 regions of three adult male brains. We identified 188 cell types and characterized their molecular signatures. Integrative analyses revealed concordant changes in DNA methylation, chromatin accessibility, chromatin organization, and gene expression across cell types, cortical areas, and basal ganglia structures. We further developed single-cell methylation barcodes that reliably predict brain cell types using the methylation status of select genomic sites. This multimodal epigenomic brain cell atlas provides new insights into the complexity of cell-type-specific gene regulation in adult human brains.

Exploring adenovirus in water environments: a systematic review and meta-analysis.

Adenoviruses (AdVs) have a significant impact in both medical and environmental contexts. The objective of this study was to investigate the prevalence of AdV in different water types, such as untreated and treated wastewater, surface water, groundwater, drinking water, and other water matrices. A total of 239 articles were included in this meta-analysis. Adenoviruses were detected in various waters worldwide. The overall prevalence in water was found to be 59.2%, with the highest prevalence in untreated wastewater (83.1%) and treated wastewater (75.3%), followed by "other water matrices" (53.4%), surface water (49.5%) drinking water (22.7%), and groundwater (18.5%). Most of the studies did not assess the viability of the viruses, leading to weak links between water contamination and risk. Both human and animal AdV were found in water environments. The findings suggest that water, including drinking water, could be a significant route of AdV transmission in both developed and developing economies.

Loss of Endothelial Glycocalyx During Normothermic Machine Perfusion of Porcine Kidneys Irrespective of Pressure and Hematocrit.

Normothermic machine perfusion (NMP) is a promising modality for marginal donor kidneys. However, little is known about the effects of NMP on causing endothelial glycocalyx (eGC) injury. This study aims to evaluate the effects of NMP on eGC injury in marginal donor kidneys and whether this is affected by perfusion pressures and hematocrits. Methods: Porcine slaughterhouse kidneys (n = 6/group) underwent 35 min of warm ischemia. Thereafter, the kidneys were preserved with oxygenated hypothermic machine perfusion for 3 h. Subsequently, 4 h of NMP was applied using pressure-controlled perfusion with an autologous blood-based solution containing either 12%, 24%, or 36% hematocrit. Pressures of 55, 75, and 95 mm Hg were applied in the 24% group. Perfusate, urine, and biopsy samples were collected to determine both injury and functional parameters. Results: During NMP, hyaluronan levels in the perfusate increased significantly (P < 0.0001). In addition, the positivity of glyco-stained glycocalyx decreased significantly over time, both in the glomeruli (P = 0.024) and peritubular capillaries (P = 0.003). The number of endothelial cells did not change during NMP (P = 0.157), whereas glomerular endothelial expression of vascular endothelial growth factor receptor-2 decreased significantly (P < 0.001). Microthrombi formation was significantly increased after NMP. The use of different pressures and hematocrits did not affect functional parameters during perfusion. Conclusions: NMP is accompanied with eGC and vascular endothelial growth factor receptor-2 loss, without significant loss of endothelial cells. eGC loss was not affected by the different pressures and hematocrits used. It remains unclear whether endothelial injury during NMP has harmful consequences for the transplanted kidney.

Integrative species delimitation in the common ophiuroid Ophiothrix angulata (Echinodermata: Ophiuroidea): insights from COI, ITS2, arm coloration, and geometric morphometrics.

Ophiothrix angulata (Say, 1825) is one of the most common and well-known ophiuroids in the Western Atlantic, with a wide geographic and bathymetric range. The taxonomy of this species has been controversial for a century because of its high morphological variability. Here we integrate information from DNA sequence data, color patterns, and geometric morphometrics to assess species delimitation and geographic differentiation in O. angulata. We found three deeply divergent mtDNA-COI clades (K2P 17.0-27.9%). ITS2 nuclear gene and geometric morphometrics of dorsal and ventral arm plates differentiate one of these lineages, as do integrative species delineation analyses, making this a confirmed candidate species.

Single-cell DNA Methylome and 3D Multi-omic Atlas of the Adult Mouse Brain.

Cytosine DNA methylation is essential in brain development and has been implicated in various neurological disorders. A comprehensive understanding of DNA methylation diversity across the entire brain in the context of the brain's 3D spatial organization is essential for building a complete molecular atlas of brain cell types and understanding their gene regulatory landscapes. To this end, we employed optimized single-nucleus methylome (snmC-seq3) and multi-omic (snm3C-seq1) sequencing technologies to generate 301,626 methylomes and 176,003 chromatin conformation/methylome joint profiles from 117 dissected regions throughout the adult mouse brain. Using iterative clustering and integrating with companion whole-brain transcriptome and chromatin accessibility datasets, we constructed a methylation-based cell type taxonomy that contains 4,673 cell groups and 261 cross-modality-annotated subclasses. We identified millions of differentially methylated regions (DMRs) across the genome, representing potential gene regulation elements. Notably, we observed spatial cytosine methylation patterns on both genes and regulatory elements in cell types within and across brain regions. Brain-wide multiplexed error-robust fluorescence in situ hybridization (MERFISH2) data validated the association of this spatial epigenetic diversity with transcription and allowed the mapping of the DNA methylation and topology information into anatomical structures more precisely than our dissections. Furthermore, multi-scale chromatin conformation diversities occur in important neuronal genes, highly associated with DNA methylation and transcription changes. Brain-wide cell type comparison allowed us to build a regulatory model for each gene, linking transcription factors, DMRs, chromatin contacts, and downstream genes to establish regulatory networks. Finally, intragenic DNA methylation and chromatin conformation patterns predicted alternative gene isoform expression observed in a companion whole-brain SMART-seq3 dataset. Our study establishes the first brain-wide, single-cell resolution DNA methylome and 3D multi-omic atlas, providing an unparalleled resource for comprehending the mouse brain's cellular-spatial and regulatory genome diversity.

Pneumococcal vaccination status among cirrhotic patients in Italy: a neglected topic.

To date, few reports have evaluated the pneumococcal vaccination status in cirrhotic patients. No data are available for European countries. We have explored this topic and the potential independent predictors motivating lack of vaccination in Italy. Between January 1st and June 30th 2022, 1419 cirrhotic patients of any etiology were consecutively enrolled in an observational, prospective study at 8 referral centers in Italy. Adjusted odds ratios (ORs) for the association with lack of vaccination were evaluated by multiple logistic regression analysis. Overall vaccine coverage was 17.9% (8.9% in patients < 65 years of age and 27.1% in those aged ≥ 65 years; p < 0.001). Among the 1165 unvaccinated patients, 1068 (91.7%) reported lack of information regarding vaccination as the reason for not having undergone vaccination. Independent predictors associated with lack of vaccination were age < 65 years (OR 3.39, CI 95% 2.41-4.76) and a higher number of schooling years (OR 2.14, CI 95% 1.58-2.91); alcoholic etiology resulted only marginally associated (OR 1.91, CI 95% 1.03-3.52). These findings establish evidence on how pneumococcal vaccination status in Italy is largely suboptimal among cirrhotic patients. These results raise concern, considering the severe outcomes of pneumococcal infection in patients with chronic liver diseases.